Symmetries of the stochastic Burgers equation 
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All Lie symmetries of the Burgers equation driven by an external random force are found. Be- 
sides the generalized Galilean transformations, this equation is also invariant under the time 
reparametrizations. It is shown that the Gaussian distribution of a pumping force is not invari- 
ant under the symmetries and breaks them down leading to the nontrivial vacuum (instanton). 
Integration over the volume of the symmetry groups provides the description of fluctuations around 
the instanton and leads to an exactly solvable quantum mechanical problem. 



PACS number(s): 05.40. +j, 31.15.Kb 

Introduction. — Statistical theory of turbulence was put 
forward by Kolmogorov in 1941 [0 and has since been 
developed intensively. The cornerstone of although sim- 
ple but surprisingly robust Kolmogorov's dimensional 
analysis is the assumption that in the fully developed 
turbulence there is a range of scales where the veloc- 
ity structure functions are universal; i.e., independent 
of the cutoffs provided by the scales of energy pumping 
and dissipation. Much efforts have been made to un- 
derstand whether there are fluctuation corrections to the 
mean field scaling exponents predicted by Kolmogorov 
and whether these corrections depend on the dissipation 
or pumping scale [0. Nevertheless, the problem is still 
far from being solved. 

Recently, the one-dimensional Burgers equation 



u t + uu x - vu xx = f(x, t) 



(1) 



driven by Gaussian random force f(x,t) with the zero 
mean and covariance 



(f(x,t)f(x',t'))= K (x-x')5(t-t'), 



(2) 



has been in the focus of quite a number of studies. The 
reason is that this equation is the simplest one that re- 
sembles the analytic structure of the Navier-Stokes equa- 
tion and, at least formally, is within the scope of applica- 
bility of the Kolmogorov's arguments. Extensive numeri- 
cal simulations of this equation || , although reproducing 
Kolmogorov's scaling exponents for the energy spectrum 
and for the two-point velocity correlation function, reveal 
strong intermittency for high-order moments of velocity 
differences. 

By adopting the hypothesis of existence of the opera- 
tor product expansion in the limit of small but nonzero 
viscosity, v — > +0, Polyakov W reduced the problem of 
calculation of high-order moments to an exactly solvable 
quantum mechanics and qualitatively explained the re- 
sults of numerical simulations. 

Another approach was based on the equivalence of the 
stochastic Burgers equation, Eqs.(|l|,^), with the Martin- 
Siggia-Rose field theory. Employing the saddle-point ap- 
proximation in the corresponding path integral, Eq. ( 14 ) , 



Gurarie and Migdal Q found the same generating func- 
tional of the velocity correlation functions as was pre- 
dicted by Polyakov. The classical solution that provides 
minimum of the action was named the " instanton" . 

The purpose of this Letter is to show that the Burg- 
ers equation driven by an external random force, Eq. (Q), 
possesses two infinite symmetry groups that are broken 
by the Gaussian distribution of a pumping force, Eq. (|^), 
thus resulting in the instanton solution. The exact in- 
tegration over the volumes of the symmetry group cor- 
responds to the description of fluctuations around the 
instanton and leads finally to the Polyakov's exactly solv- 
able quantum mechanics. 

Martin-Siggia-Rose formalizm. — Let us first consider 
the functional Langevin equation of the general type 



d t u + T(x, [u]) = f(x, t), 



(3) 



here f(x,t) is a functional Gaussian process with the 
covariance operator K. Local functional !F{x, [u]) com- 
pletely characterizes the dynamics of the system. 

Then, the Martin-Siggia-Rose formalizm permits one 
to obtain a formal path integral representation for the ex- 
pectation values of the observables 0[u] in the stationary 
state § [K} 



(0[u\) 



VfO[u] 5{d t u + T-f)e 



-(f,K- 1 f)/2 



(4) 



VfWf 0[u] 8{d t u + T- f)e 



where, we have introduced an auxiliary field fi to rewrite 
Gaussian path integral over the field /. Now, let us pass 
in this path integral from the integration over the force 
field Vf to the integration over the velocity field Vu. 
Using the (^-function we obtain 

(0[u\) = J VyVu 0[u) J[u) exp(-S[u,n])- (5) 

The action of the effective field theory can be written as 

S[u,fjb] = -i(jj,,d t u + J T ) + '^{fj,,Kfi). (6) 
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The functional J[u] is the Jacobian of the force-to- 
velocity transformation 



can be represented by the path integral 



J[u] = det 



'11 
Su 



= det 



d t + 



8JF_ 
Su 



(7) 



Here, the variational derivative ST / Su is determined by 
the usual equality 

F{x, [u + Su}) - T(x, [u]) = / S:F } X : ^ Su{y) dy. (8) 



Su(y) 

Formal calculation of the determinant, Eq. (fjj), leads to 
the result 



J[u] = exp ( 0(0) Tr -=- 
ou 



(9) 



involving ill-defined quantity 0(0), where 9(t) is usual 
step function. To give a meaning to this expression one 
can make use of the usual limiting procedure in which 
time interval (U,tf) is subdivided into TV equal subin- 
tervals of duration At each, and the differential relation 



inside is replaced by a difference one |10|. After taking 
the continuous limit one finds: 9(0) = 1/2. 

In the perturbative treatment of the path integral this 
Jacobian term is irrelevant because it can always be ab- 
sorbed by the proper redifinition of the Green functions. 
Nevertheless, perfoming manipulations with the path in- 
tegral we should take care about this Jacobian term. For 
stochastic Burgers equation the local functional T(x, [u]) 
can be defined as follows: 

F[x, [it]) = J dy 8(x - y) (uu y - vu yy ) . (10) 

Its variational derivative, Eq. (^), is then equal to 
8T(x, [u]) 



= S'(x-y)u-vS"(x-y), 
5u(y) 

and the trace term formally becomes 

Tr^ = fdtdx ^M 
Su J Su(x) 

= S'(0) [ dtdx it- (0) / dtdx. 



(11) 



(12) 



Here, the first term vanishes for symmetry reasons and 
the second term does not depend on u and can be ab- 
sorbed by the normalization factor in the path integral. 
Hence, for the Burgers equation the Jacobian of the 
force-to-velocity transformation is equal to some irrele- 
vant constant. 

So, the generating functional of the velocity correlation 
functions 



(F[X\) = ( exp / dx X(x)u(x,0) ) , / dx X(x) = 0, 

(13) 



J T>(iVuF[\] exp(-S[u,fi]), 



(14) 



where action of the effective field theory is equal to 
S[u, u] = — i / dtdxfi(u t + uu x — vu xx ) 



(15) 



+ — / dtdxdx' n(x,t)K,(x — x')/j.(x' ,t), 

and time integration runs over the interval (£j,t/). For 
the sake of simplicity, we will consider the time interval 
(-co, 0) although all the following does not depend on 
the particular choice. The covariance k(x) is supposed 
to be a slowly varying even function with the expansion 



k(x) = k(0) - y X 2 , 



M « 



'k(O) 

K 



I. 



(16) 



and quickly turns into zero when |x| 3> I. The interval I 
characterises the correlation length of the random force 
and we only study velocity correlation functions within 
this interval. 

Symmetries of the Burgers equation. — We start our 
consideration with a quite natural question: what is the 
most general transformation leaving invariant the Burg- 
ers equation driven by an external force, Eq. ([!]), or, in 
other words, given two functions u(x,t) and f(x,t) sat- 
isfying the Burgers equation, what is the most general 
transformation that produces another pair of functions 
satisfying the same equation? 

Making use of the methods of Lie group analysis of dif- 
ferential equations [ p~T| it is possible to prove that all such 
transformations form the symmetry group which consists 
of the two infinite subgroups: generalized Galilean trans- 
formations (G), 

t = t, (17a) 
x = x + a(t), (17b) 
u = u + a'(t), (17c) 
f = f + a"(t), (17d) 

and time reparametrizations (L), 

i=b(t), (18a) 
x = xy/Vtt), (18b) 

/=(/ + |Wi),f})/vWP. (18d) 

where a(t) is an arbitrary function; b(t) maps the time 
interval (ti,i/) onto itself and b'(t) > 0. The braces 



ri , , , b'" 3 fb" 



(19) 
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are known in the theory of complex functions as the 
" Schwarzian derivative" . 

Symmetry transformations (|lj],[l8|), although not de- 
pending explicitly on the viscosity ^, are stipulated by 
the structure of the diffusion term. In particular, these 
transformations leave also invariant, along with the Burg- 
ers equation itself, the following relation, 



(dx — udt) 2 ~ dt, 



(20) 



which coincides (up to convective term, udt) with a sim- 
ilar relation for a diffusion process. 

The symmetry groups corresponding to the transfor- 
mations (|l7|,[l8|) are generated by the infinitesimal oper- 
ators 



d ^ d , a /id 
dx du df ' 



G(a(t)) = ct— — -|- n — i- n 



(21a) 



L ^=4 + ^'4x + 1 2^ X -^l 

+ l((3'" x -3(3"f)^j, (21b) 

which, in the basis 

G r = G(f +1 l 2 ), L n = L(t n+1 ), (22) 
form Lie algebra with the commutation relations 



[G r ,G s ] = 0, 

[L n ,G 8 ] = (s-n/2) G n+S , 
[L n , L m ] = (to — n) L n + m , 



(23a) 
(23b) 
(23c) 



here n, m are integers and r, s are half-integers. 

Action transformation laws. — Up until now, we have 
not been concerned about the probability distribution 
of the pumping force. However, if the force f(x,t) is a 
Gaussian random function, then its probability measure 
is not invariant under the symmetry transformations. It 
is easy to check directly that the action ( fL5| ) which in 
fact defines such a probability measure changes under 
the generalized Galilean transformations as 



S[u, p] = S[u, fj]-i J dtdx fi(x, t) a"(t), (24) 
and under the time reparametrizations as 
S[u, p] = — i / dtdx[j,(ut + uu x — vu xx ) 



+ - I dtdxdx' b' 2 (i(x,t)n{Vb'(x - x'))fi(x',t) 



— J dtdx xfi(x,t) {b(t),t}. 



(25) 



Hence, we come to the conclusion that it is the Gaussian 
distribution of the pumping force which breaks down the 
infinite symmetry group of the Burgers equation (^) . 



The action is still invariant under the finite subgroup 
which consists of spatial translations generated by the 
infinitesimal operator G_i/2, Galilean transformations 
generated by the operator G1/2, and time translations 
with the generator L^\. According to the Noether theo- 
rem [|ll] there are three conservation laws corresponding 
to the three-parameter subgroup of variational symme- 
tries: the momentum conservation 

P = —i J dx n u x , (26) 

the conservation of the center-of-mass motion 

M = i I dx fj, + tP, (27) 



(28) 



here L is the Lagrangian S = j dtL. 

The Faddeev-Popov method. — If the action were invari- 
ant with respect to the infinite symmetry group, we could 
employ the Faddeev-Popov method to eliminate the cor- 
responding degrees of freedom p!o[ . As we will see later, 
in our case the same method separates the degrees of free- 
dom of the symmetry group transformations and leads fi- 
nally to the exactly solvable quantum mechanics for the 
separated modes. The Faddeev-Popov method consists 
of the three steps. At first, considering the field theory 



and the energy conservation 

E = i / dx fi Ut + L, 



VAexp(-S[A]), 



(29) 



with the action, S[A], invariant under the group Q of 
gauge transformations S 1 ^ 9 ] = S[A], we define the equa- 
tion f(A) = which fixes uniquely the gauge degrees 
of freedom; i.e., the equation f(A 9 ) — should provide 
the unique solution for the group element g for an arbi- 
trary chosen field configuration A. Then, we define the 
functional Jacobian Jg [A] by the condition 



Jg[A] / Vg8(f(A°)) = l ) 



(30) 



where Jg[A] = det \\Sf(A 9 )/Sg\\, and integration runs 
over the volume of the group Q. It is clear that the func- 
tional Jg [A] is invariant by construction under the gauge 
transformations, Jg[^4 9 ] = So, we can calculate 

the Jacobian only for the identity element of the group. 
Finally, if we insert the identity (30) into the functional 
integral ( ^9| ) and shift the variables A 9 — > A, we sepa- 
rate the degrees of freedom corresponding to the group 
volume 

Z = JvA J g [A] 6(f(A)) exp(-S[A]) J Vg. (31) 
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Now, we can apply the same method to integrate over 
the volume of the groups of generalized Galilean trans- 
formations and time reparametrizations. 

Integration over the volume of the group G. — The 
gauge fixing equation in this case is: m(0, t) = 0, and 
we consider the identity 



J G [u] / Va(t) 5 {a - u(a, t)) = 1, 



(32) 



which is the definition of the Jacobian Jg [u] ; integration 
runs over the volume of the group G. The explicit form of 
the Jacobian, Jq, is given by the functional determinant 



Jg[u] = dct 



di 



Su(a(t),t) 



5a{t) 



(33) 



For the identity element of the group, a(t) = 0, this ex- 
pression takes the form 



JbN = dfit||at-u a (0 ) t)|| 



(34) 



To preserve the causality of the Langevin equation after 
the generalized Galilean transformation we should prop- 
agate operator, whose determinant we are calculating, 
backward in time. Standard calculation [fl0|| then gives 
the result 



Jq[u] = exp (1 - 6(0)) / dt u x (0,t) 



(35) 



involving the same ill-defined quantity 8(0) that appeared 
in the derivation of the Martin-Siggia-Rose field theory, 
Eq. (||). As we will see later the choice, 8(0) = 1/2, 
is self-consistent in the sense that it ensures finally the 
existence of the steady state. Actually, the value of 8(0) 
plays exactly the same role as the B-anomaly term in the 
Polyakov operator product expansion. 

Now, we insert the identity ( |32| ) into the functional 
integral (|l3|). Then, after the properly defined general- 
ized Galilean transformation ( |l7| ) , we turn the argument 
of the 5-function into the gauge fixing term; i.e., turn it 
into the equation of the surface in the functional space 
that intersects the orbits of the symmetry group only 
once. Finally, we get for the generating functional of the 



velocity correlation functions (13) 



(F[X}) = Jvtijva exp(-S G [a,fj]) 

Vu 5(u(0)) J G F[X] exp(-S[u,/j]), (36) 



here we used the condition, J dx X(x) = 0, which ensures 
the invariance of the functional F[X] under the general- 
ized Galilean transformations; 



Sg[o>, A*] = — * J dtdx n(x,t) a"(t) 



(37) 



is the effective action of the corresponding modes. In 
fact, all these modes can be integrated out, leading to 
the additional constraint 



(38) 



Va(t) exp(— Sc[a, jj]) = S I d t I dxji 



From this it follows that the zero moment of the auxiliary 
field fi is just a linear function of time 



J dx/j, = i(Pt - M), 



(39) 



where P and M are the integrals of motion (|26|,[27|). In 
the frame of reference moving together with the center of 
mass of the Burgers fluid both of the integrals are equal 
to zero and itq = 0. 

The integration over the volume of the group of gener- 
alized Galilean transformations has first been proposed 
in Refs. PJl^] to justify the stability of the saddle-point 
approximation in the path integral. 

Integration over the volume of the group L. — The 
gauge fixing equation is: u(xQ,t) = 0, where xq ^ 0, 
and we consider the identity 



J L [u] J Vb(t) 5 (^y - u(x Vv,b) VV 

The explicit expression for the Jacobian Jl is 
5 



= 1, (40) 



Jl[u] = dct 



which in the new variables 

can be rewritten as 

S {x a' - u(x e°, b)e- a } e 2a 



(41) 



(42) 



Jl[u] = dot 



5a 





5a 


det 




~Sb 



(43) 



Again, due to invariance of the Jacobian under the group 
of time reparametrizations it can be computed only for 
the identity element of the group 



J L [u] = det \\d t - u x (x Q ,t)\\ 



(44) 



where the unessential constant coefficient that does not 
depend on u is omitted and constraint, u(xo,t) = 0, is 
used. Calculating the determinant along the lines indi- 
cated above we get 



J L [u] = exp Q J 



dt u x (x ,t) 



(45) 



Let us now insert the identity (^) into the functional 
integral (|3q). Then, after a properly defined time 
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reparametrization, the generating functional of the ve- 
locity correlation functions is converted into 



Vu 6{u(0)) S(u(x )) J G J L F[X] exp (S ) 



J Vb exp j-S L [M] + j dx x\(x) \ (4(i) 

where ir{t) — J dx xfi(x,t) and 



So[u, n] = -i / dtdxfi(u t + uu x - vu xx ), 



(47) 



is the action of the modes remaining. Due to the special 
choice of the noise covariance, Eq.([l6|), it does not depend 
on the pumping force. 

In the limit of the small but nonzero viscosity, v — > +0, 
we can separate also the integration over tt. The desired 
effective action of the separated modes b(t) and tr(t), 



S L [b,ir] 



2 

K0 



dt w(t) {b{t),t} - 
dt b'(tf ir(t) 2 , 



dt 



b"(t) 
2b'(i) 



(48) 



follows then from the action transformation law under 
the time reparametrizations, Eq. (|25|), and the analogous 
transformation law for the Jacobian 



If b" 

J G [u] = J G [u] exp ( - / dt — 



(49) 



Separating the modes, we obtain another ambiguous 
parameter. Namely, we can choose an arbitrary normal- 
ization factor, Zq, for the partition function of the modes 
of the symmetry group transformations and the inverse 
one, 1/Zq , for the partition function of the modes remain- 
ing. This parameter plays exactly the same role as the 
A-anomaly term in Polyakov theory. Again, the choice 
Zq = 1 is self-consistent and it ensures the existence of 
the steady state (see also [[u|). 

It should be stressed at this point that the gauge fix- 
ing term 6(u(xq)) works properly only if the degrees of 
freedom associated with the field a(t) have already been 
integrated out; i.e., we cannot change the order and in- 
tegrate first over the volume of the group L. 

Equivalent Quantum Mechanics. — The action ( f48| ) is 
equivalent to the Polyakov's exactly solvable quantum 
mechanics. To make it obvious, let us pass from the 
variables t, b(t) and ir(t) to the new ones 



'' = '""■ </ !/ ' i = S^ P(b) = b'(t)n(t), 



2b'{tf 



and finally get 



S L [q,p] = J db i-ip (q 



2\ ^ Kq 2 

q) -2 q + Y p 



(50) 



(51) 



where the additional term, —q, comes from the Jacobian 
T>itT>b/T>pT>q. This action corresponds to the quantum 
mechanics (in imaginary time) with the Hamiltonian |l(J 



2 V K o 



g 

2k 



3q 
2 ' 



(52) 



Following Polyakov ||] , we can find the zero energy eigen- 
function of this Hamiltonian and finally calculate the gen- 
erating functional of the velocity correlation functions, 



(F[A]> = exp • 



dx x\(x) 



3/2" 



(53) 



In conclusion I would like to mention that the method 
proposed in this Letter seems to have a wider applica- 
bility than just to the Burgers equation considered. Any 
stochastic equation of the Langevin type that possesses 
infinite symmetry group can be treated this way. 
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